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Abstract 

We calculate the baryon flow (P x /A(y)) in the energy range from 
.25 to <2.5AGeV in a relativistic transport model for Ni + Ni and 
Au + Au collisions employing various models for the baryon self ener- 
gies. We find that to describe the flow data of the FOPI Collaboration 
the strength of the vector potential has to be reduced at high relative 
momentum or at high density such that the Schrodinger- equivalent 
potential at normal nuclear density decreases above lGeV relative 
kinetic energy and approaches zero above 2GeV. 
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1 Introduction 



The nuclear equation of state (EOS) at high density (p > 3po) is still an 
unsettled issue though many experimental efforts have been made in the 
last couple of years to address this question in a more systematic way ]!]. |]. 
[|. [3j]. Experimentally, the baryon sidewards flow and subthreshold particle 
production are the most promising observables. Recently |3|, |5|, the flow 
has been measured for heavy- ion collisions at SIS energies (< 2AGeV) for 
Ni + Ni and Au + Au systems, which provides further constraints on the 
hadronic models and the EOS employed. 

From the theoretical point of view relativistic transport models || [7|. || 
have been used to describe heavy-ion collisions at energies from the SIS at 
GSI to the AGS at BNL and the SPS regime at CERN (< 200 AGeV). 
There are two main ingredients in the transport model: the mean fields (i.e. 
scalar and vector self energies) and in-medium baryon-baryon and meson- 
baryon cross sections. The scalar and vector mean fields are usually derived 
from some hadronic Lagrangian (which gives a well defined EOS) while the 
baryon-baryon and meson-baryon cross sections are taken from experimental 
data for the related processes in free space. In-medium modifications of the 
inelastic channels, furthermore, are constrained by the experimental data on 
the particle rapidity distributions which control the amount of 'stopping'. 

In principle, the baryon selfenergies should be determined by a Dirac- 
Brueckner approach including all relevant hadronic couplings. However, cal- 
culations for configuration dependent phase-space distributions are not very 
reliable yet and limited in density as well as in momentum |5[ |H], TT|, [12| . 

As mentioned before, the scalar and vector self energies for nucleons with 
their explicit momentum and density dependence are the key quantities that 
determine the nuclear EOS. In this work we will perform a systematic study 
of Ni + Ni and Au + Au collisions in order to extract further information 
on these quantities in comparison to the recent experimental data on the 
collective flow of baryons. 

Our work is organized as follows: In Section 2 we will briefly describe the 
relativistic transport approach employed as well as the known constraints on 
the momentum dependence of the scalar and vector self energies. Section 3 is 
devoted to a systematic comparison of the calculated flow - employing various 
Lagrangian models - to the experimental data while Section 4 concludes with 
a summary and discussion of open problems. 
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2 The transport model 



In this work we perform the theoretical analysis along the line of a relativistic 
transport approach which is based on a coupled set of covariant transport 
equations for the phase-space distributions fh(x,p) of a hadron h J7|, [13[] , i.e. 

{(n„ - n„flgt£ - M*^c/f) ^ + (ru^ + m*^) 9;} /,(x, P ) 
= £ / rf2rf3cM [c f G'] 12 ^34 #(n + n 2 - n 3 - n 4 ) 

x { A 3 Ps) A 4 Pdfh(x, p)fh 2 (x,p 2 ) 

- fh(x,p)f h2 (x,p 2 )fh 3 (x,P3)f h4 (x,p 4 )} . (1) 

In Eq. ([I]) Uft(x,p) and U£(x,p) denote the real part of the scalar and vector 
hadron self energies, respectively, while 12^34^ (IT + Yl 2 — II3 — II4) is 

the 'transition rate' for the process 1 + 2 — > 3 + 4. Though in quantum many- 
body systems the transition rate is partly off-shell - as indicated by the index 
r of the 5-function - we use the semi-classical on-shell limit r — > since this 
approximation is found to describe reasonably well hadronic spectra in a wide 
dynamical regime. The hadron quasi-particle properties in (|l|) are defined 
via the mass- shell constraint 5(11^ - Mf) g3| with effective masses and 
momenta given by 

M*(x,p) = M h + U%(x,p) 

U»(x,p) = p^-U^x,p) , (2) 

while the phase-space factors 

fh(x,p) = 1 - f h (x,p) (3) 

account for fermion Pauli-blocking. The transport approach ([!]) is fully spec- 
ified by Ufr(x,p) and Uj*(x,p) (// = 0, 1, 2, 3), which determine the mean-field 
propagation of the hadrons, and by the transition rates G'G8^(. . .) in the 
collision term that describe the scattering and hadron production/absorption 
rates. For the latter we employ free cross sections as in Ref. |14] that are 
parameterized in line with corresponding experimental data. In the relativis- 
tic transport approach we explicitly propagate nucleons and A's with their 
isospin degrees of freedom. For more details we refer the reader to Ref. jlSj . 
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Before going over to a discussion of the scalar and vector self energies we 
start with the cascade limit II s = U u = and compare the flow F defined 
by 

F = ^(PM)y=o , (4) 

where y is the nucleon rapidity in the center-of-mass system normalized by 
the projectile rapidity, i.e. y = y cm /y P roj- Fig. 1 shows the calculated flow 
F (divided by A{ /3 + A 1 ^) function of beam energy for Ni + Ni and 
Au + Au systems. We first observe that in the cascade mode we obtain a 

A 1/3 A 1/3 

scaling of the flow F with the system size expressed by A\ + A 2 for both 



the BUU []Tb| and RBUU model, which will be solely used later on. It can 
be observed that the two models agree very well below lAGeV, whereas 
they differ in the energy range above. This difference comes mainly from 
the fact that the BUU includes all higher resonances up to a mass of 2GeV 
whereas the RBUU treats only the A(1232) explicitely. These high mass 
resonances, though only present in a small percentage of all particles in the 
high density phase, contribute significantly to the flow at high energies. In 
order to correct the RBUU results for this effect when including potentials 
we add the difference between BUU and RBUU in the cascade mode to the 
RBUU results. This correction of course neglects higher order effects caused 
by the smaller velocity of a heavier resonance in the surrounding nuclear 
matter; nevertheless, we expect these corrections to be small. 

The flow generated without any mean fields is due to baryon-baryon col- 
lisions and underestimates the experimental data for Ni + Ni at all energies 
considerably as shown in Fig. 1. The difference between the cascade results 
and the data thus is caused by the baryon self energies which we try to 
simulate in the following. 

In order to compare our flow values with other theoretical groups, we 
display the quantity 

1 y T dN 



(P x /N) dir = - J dy ( Px /N)(y) ^- sgn(y) (0) 



Vcm 



versus beam energy in Fig. 2 for Au + Au systems. The UrQMD calculations 
of the Frankfurt group predict the flow for Au + Au collisions to increase 
with bombarding energy at least up to 4AGeV |T7] while hydrodynamical 



calculations including only hadronic matter predict a decrease of the baryon 
flow above 4 — 5AGeV JTS] for the heavy system Au + Au. Furthermore, 
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hydrodynamical calculations including a transition to a quark-gluon-plasma 
(QGP) phase predict a vanishing baryon flow at 4 — 5AGeV [18[]. In Fig. 2 
we also show the results of RBUU (hard), BUU (hard) and UrQMD (hard) 
EOS calculations as well as our BUU cascade calculation in comparison with 
a UrQMD cascade calculation. The assignment 'hard' here corresponds to 
a nuclear incompressibility K « 380 MeV. We see, first, that all transport 
models agree quite well when using a 'normal' hard equation of state as well 
as in the pure cascade mode. All models predict an increasing baryon flow 
with bombarding energy. Second, the calculation with a hydrodynamical 
model shows an even stronger flow at the energies considered here. Third, 
the RBUU results with a special momentum- dependence of the mean fields 
(NL3*, see below) show a decreasing flow above lAGeV. This specific model 
will be discussed in more detail later. 

In Fig. 3 we plot the flow F (divided by A 1 / 3 + A 1 / 3 ) clS db function of the 
beam energy for both hard and soft EOS without explicit momentum depen- 
dence of the potentials. The soft EOS is taken from the work of Furnstahl et 
al.fl9|j with a nuclear incompressibility K ~ 194 MeV while the hard EOS is 
the NL3 parameter set from [|14]] with a nuclear incompressibility K = 380 
MeV. We notice from this figure that the baryon flow in the RBUU approach 
does not sensitively depend on the nuclear incompressibility. Therefore, in 
the following we focus on the momentum dependence of these potentials only. 

In another RBUU approach || |2(| calculations have been done recently 
for heavy-ion collisions by employing density dependent scalar and vector 
potentials self-consistently as well as the momentum dependence of these 
potentials while taking care of chiral symmetry constraints ||19|| . In these 
approaches the flow has been calculated also for Ni+Ni and Au+Au systems 
and been found to overestimate the experimental data at 1.5 and 2AGeV 
for Ni + Ni considerably. Thus also the latter transport calculations yield 
an increasing (or at least constant) flow with bombarding energy, whereas 
the experimental data indicate a decrease for Ni + Ni above lAGeV. We 
will argue that this decrease of the flow F puts stringent constraints on the 
momentum- and density dependence of the mean fields. 

The model inputs for the mean fields are related to the nuclear incom- 
pressibility K at density po as well as to the momentum dependence of the 
mean fields as first pointed out by Gale et al. |^1[ and incorporated later 



on also in relativistic transport models by several authors || In the 

RBUU approach - due to covariance - the scalar and vector mean fields have 
to be explicitly momentum dependent |13| in order to describe properly the 
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Schrodinger-equivalent optical potential denned by 



U sep (E km ) = U S + U + ^{U 2 - U 2 ) + ^E kin (6) 

as a function of the nucleon kinetic energy E^m with respect to the nuclear 
matter rest frame. However, above Ekm = lGeV the Schrodinger-equivalent 
optical potential is not well known experimentally, such that the flow data 
from the FOPI Collaboration could provide further constraints also on this 
quantity. 

In this work we use a similar Lagrangian density as proposed by Walecka 
|| for the description of nuclear matter, which has been used in the rel- 
ativistic BUU model before WM. This Lagrangian contains nonlinear self- 
interactions of the scalar field U(a) = ^m 2 a 2 + \Ba z + \Ca A where the 
parameters m CT , B, C are calculated by fitting the saturation density, binding 
energy, effective nucleon mass as well as the compression modulus at nuclear 
matter density (cf. NL3 parameters set from Table 1 in Ref. [T3[]). An addi- 
tional coupling between the vector and scalar fields in the Lagrangian leads 
to a relatively soft EOS [TP|, |23]. In the present calculation we do not consider 



this effect, first, because our main object is to concentrate on the momentum 
dependence on the vector and scalar fields in comparison to the flow data, 
second, because the influence of the momentum dependence is much stronger 
than that of the density or incompressibility as pointed out before. 

The energy density in mean-field theory (for nuclear matter) in these 



models can be written as 14 



e(m*,n b ) = g v V n b - \m 2 V% + ^(m- m*) 2 + -^(m - m*) 3 

C , f k f d 3 p n— r . , 

+ 4^l(™-™) +7/ o (^vV + ™*), (7) 

where m* = m — g s a is the effective nucleon mass, n b is the baryon density 
and the spin and isospin degeneracy is 7 = 4. o and Vq are the scalar and 
vector fields with mass m a and m v , which couple to nucleons with coupling 
constants g s and g v , respectively. The quantities B and C are constant 
parameters and p is the nucleon momentum which has to be integrated up 
to the fermi momentum kf. In this model the vector and scalar potentials 
are density dependent, however, the vector potential increases only linearly 



with density. Recently [20] this Lagrangian has been extended - maintaining 
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chiral symmetry constraints - to include a nonlinear dependence of the vector 
potential with density, too. We found, however, that this type of Lagrangian 
density with momentum dependent fields underestimates the flow data at all 
beam energies for Ni + Ni as well as Au + Au systems. The calculated flow 
values follow the cascade results as the scalar and vector mean fields cancel 
each other approximately for such type of EOS (cf. Fig. 1). 

In our present calculation, we use the energy density Eq.(|7|) for calculating 
the scalar and vector potentials as a function of density. Momentum depen- 
dent potentials are obtained by fitting the Schrodinger equivalent potential 
(IB) according to Dirac phenomenology for intermediate energy proton-nucleus 
scattering ||22| . We use the approach from Ref. |7j to take into account this 
momentum dependence by introducing additional scalar and vector cutoffs 
A s , A„. The scalar and vector form factors at the vertices are 

^-(p-(p)) 2 anH K-(p-(p)) 2 (R) 

A 2 + {p _ {p)) 2 A 2 + {p _ {p)) 2 > I ) 

respectively, where (p — (p)) accounts for the difference of the one-particle 
momentum to the average momentum of the surrounding nuclear matter. 
The values of A s and A„ vary from 0.95 to 1.05 GeV and 0.9 to 1.0 GeV to 
get a good fit to the data. This momentum dependence is not computed self- 
consistently on the mean field level since it leads only to a small change in the 
original parameters of the model as well as in the fitting of the Schrodinger 
equivalent potential. So this approximation does not effect much our flow 
results, because for nuclear matter the energy scale involved is much smaller 
(kp <C A s , A„) than in the initial stage of high-energy heavy-ion collisions. 

The Schrodinger equivalent potential (§) is shown in Fig. 4 as a function 
of the nucleon kinetic energy with respect to the nuclear matter at rest; 



also plotted are the data from Hama et al. p2[ . The solid line in Fig. 4 

-jp-(p)? 
<p-(p)) 2 



is for the special momentum dependence including form factors 



and to ; 1 P m?2 ■ The dashed line is obtained using the scalar and vector form 

K + YP-KP)) to 
A 2 A 2 

factors A 2_^ p2 and A2 ^ p2 with A s = 0.95 GeV and A^ = 0.9 GeV, respectively. 



3 Comparison to experimental data 

We use the same parameter sets as for the Schrodinger equivalent potential 
in our flow calculations. The calculations are performed for the impact pa- 
rameter b = 4fm for Ni + Ni and b = 6fm for Au + Au systems, since for 
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these impact parameters we get the maximum flow, which corresponds to 
the multiplicity bins M3 and M4 as defined by the Plastic Ball collaboration 
PHI . We have calculated the flow, i.e. the slope parameter (Q) by fitting a 
straight line from -0.3< y < 0.3 for Ni + Ni and Au + Au systems at all 
energies. Higher order terms (e.g. fitting a polynomial of 3rd grade) didn't 
change the results systematically. 

In Fig. 5 and 6 the flow F (divided by (A^ 3 + A 1 ^ 3 ) as in Fig. 1) is 
displayed in comparison with the data from Refs. 0, || for different sys- 
tems. Fig. 5 shows the flow for Ni + Ni in the energy range up to 2.5AGeV 
and Fig. 6 for Au + Au, respectively. In both Figs. 5 and 6 the solid line 
(NL3) is obtained without explicit momentum dependence of the self ener- 
gies, whereas the dashed line corresponds to the momentum dependent scalar 
and vector potentials and the dashed-dotted line (NL3*) corresponds to the 
special momentum dependence shown in Fig. 4. 

From Fig. 5 we observe that the dashed-dotted line (NL3*) is in good 
agreement with the flow data for Ni + Ni, whereas for Au + Au (Fig. 6) 
the dashed-dotted line (NL3*) is slightly below the flow data from []3], |5| but 
closer to the plastic ball data (open triangles). 

The important point to be noted here is that the flow rises up to lAGeV 
and decreases above lAGeV for Ni + Ni, whereas it saturates above lAGeV 
for Au + Au. The physics behind is that the repulsive force due to the vector 
mean fields must decrease considerably at high beam energy such that the 
Lorentz force on the particles vanishes in the initial phase of the collision. In 
subsequent collisions, which are more important in the Au+Au case due to its 
size, the kinetic energy of the particles moving relative to the local rest frame 
is then in a range where the Schrodinger equivalent potential is determined 
by |p2[ . We thus conclude that to explain the flow data up to 2AGeV one 
has to reduce the vector mean field considerably. In other words, there is 
only a weak repulsive force at high relative momenta and high densities. 

Another interesting point is that the flow in our calculation for Ni + Ni 
decreases earlier with beam energy than for the Au + Au system. This is due 
to fact that the flow is governed by the average transverse pressure generated 
due to the number of nucleons. For these reasons the flow for Au + Au first 
saturates and then decreases at a higher beam energy (> 2.5AGeV). This 
implies that the observed A 1 / 3 scaling of the flow up to lAGeV does no 
longer hold for higher energies. Here it would be very interesting to have - 
besides the Ni + Ni system recently measured at FOPI - systematic data 
for higher mass systems and beam energies above lAGeV. For Au + Au 
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first preliminary results have been presented by the EOS Collaboration [25 
indicating a gradual decrease of the flow F from 2 - 8AGeV. 

Finally, we show in Fig. 7 the EOS as well as the scalar and vector po- 
tential energy associated with the special momentum dependence (NL3*). 
The upper part shows the energy per nucleon in comparison to the standard 
NL3 parametrization |TJJ] (dashed line). For NL3* (solid line) the approxi- 
mate equation of state has been derived by fitting the vector potential with 
an appropriate polynomial in the baryon density. The lower part shows the 
corresponding vector and scalar potential energy as a function of the baryon 
density. The vector part for NL3* is substantially lower at high baryon 
density as compared to the NL3 parameter set as well as to the original 
o — cu-model ||. The vector potential at 3p is about 460, 653 and lQ2QMeV 
for NL3*, NL3 and the o — uj- model, respectively. 



4 Summary 

In this work we have calculated the baryon flow in the energy range up to 
2.bAGeV in a relativistic transport model for Ni+Ni and Au+Au. We found 
that in order to properly describe the flow data of the FOPI Collaboration 
at high beam energies, the strength of the vector potential has to be reduced 
considerably in the RBUU model at high relative momenta and/or densities. 
Otherwise, too much flow is generated in the early stages of the reaction and 
cannot be reduced at later phases where the Schrodinger equivalent potential 
is experimentally known and constrained. 

This assumed decrease of the vector and scalar potentials destroys, how- 
ever, the A 1 / 3 scaling of the flow - observed in heavy ion collisions below 
lAGeV - at higher energies. Therefore it would be interesting to have sys- 
tematic studies for different mass systems for higher bombarding energies. 
From these studies it might even be possible to extract the dependence of 
the potentials on density and momentum separately. 

One shortcoming of the transport model used here is the restriction to 
binary nucleon-nucleon or meson-nucleon scattering. Especially in meson- 
nucleon reactions at 2AGeV bombarding energy the baryon excitations be- 
come very high (a/s pa 2.5GeV) such that multi-particle final states can occur 
which may lead to a decrease of directed flow. Despite of these uncertainties 
above ~ 2AGeV, the decrease of the Schrodinger equivalent potential above 
lAGeV should have clearly been demonstrated. 
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Fig.l The flow F{y) versus beam energy per nucleon for Ni + Ni and impact 
parameter b = 4 fm and Au + Au systems with impact parameter b = 6 fm in the 
cascade mode. The solid line corresponds to the RBUU results for Au + Au, the 
dashed line to the RBUU results for Ni + Ni, whereas the dotted and dashed- 
dotted lines correspond to the BUU results for Au + Au and Ni + Ni systems, 
respectively. The experimental data for Ni + Ni are from Ref.Q. 



Fig. 2 (P x j 'N) dlT '(y) values versus beam energy per nucleon for Au+Au at b = 4fm. 
A hydro dynamical calculation (for b = 3 fm) is shown versus UrQMD (hard EOS), 
BUU (hard momentum dependent EOS) and RBUU with momentum dependence 
(NL3) according to Q. Also shown are the UrQMD cascade calculations from 
Ref.|l8| and the BUU cascade calculations pf, which are found to agree very 
well. The dashed line results from a RBUU calculation with special momentum 
dependence (NL3*). 

Fig. 3 The flow F{y) versus beam energy per nucleon for Ni + Ni at b = 4 fm. 
The dashed line results from the soft EOS from Ref. [|HJ while the solid line is 
obtained for a hard EOS [14|. The data points are from the FOPI Collaboration 



Fig. 4 The Schrodinger equivalent potential (^) as a function of the nucleon kinetic 
energy Efo n . The solid curve is for the special momentum dependence NL3* dis- 
cussed in the text and the dashed curve for the momentum dependent parameter 



set NL3 M (see text). The data points are from Hama et al. [22]. 



Fig. 5 The flow F(y) versus the beam energy per nucleon for Ni + Ni at b = 4 
fm. The solid line results for the parameter set NL3, the dashed line includes 
an explicit momentum dependence of NL3 and the dashed-dotted line is for the 
special momentum dependence NL3*, where the vector potential decreases at high 
energies (c.f. Fig. 4). The data points are from the FOPI Collaboration Q. 

Fig. 6 Same as Fig. 5 for Au + Au at b = 6 fm. The data points are from the Plastic 
Ball, FOPI and EOS Collaborations, Ref.g. 

Fig. 7 Upper part: The solid line shows the equation of state for the parameter set 
NL3* in comparison to NL3 (dashed line). Lower part: Vector (upper solid line) 
and scalar (lower solid line) potential energy per nucleon for NL3* in comparison 
to NL3 (dashed line). 
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